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ABSTRACT 


This thesis reviews four of the existing matrix methods for finding 
the shortest path in a network, including the little known matrix method 
by Floyd. Floyd's method is then extended to determine all best paths. 
After a brief review of the nth best path problem, Floyd's method is at 
again extended to determine the nth best path. Finally, the nth best . 
path problem is investigated to determine its applicability to the travel- 


ing salesman problem. 
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1. Introduction 

Four matrix methods to determine the shortest path in a network 
wit be presented. After an explanation of each method there will follow 
a simple example to illustrate the procedure. | 

First, the method of matrix manipulation [3] will be set forth. 
Next, an improvement on this, the Cascade algorithm [7], will be pre- 
sented. The two algorithms that appear to be the most efficient are 
next. Algorithm 97 was published by Floyd [8] in 1962 and went unnoticed 
for some time. In 1965 Murchland [15] published what he called “a new 
method," but which in reality was the same as Floyd's. A different 
approach but equally as efficient as Floyd's was put forth in 1966 by 
Dantzig [6] and is an inductive procedure. 

In many practical problems a shortest route is not sufficient. All 
alternate shortest routes or the n best routes are needed. Thus a 
simple extension to Floyd's algorittim to determine all alternate routes 
will be presented. Following this will be a brief review of the nth best 
path problem and a second extension to Algorithm 97. This extension will 
determine the n best paths from all nodes i to all nodes j. 

A logical extension of this method of determining all n best paths 
is to the traveling salesman problem. Although this problem is not solved 


here, a possible approach is set forth. 


2. Notation 
A graph or network, denoted G = (X,U), is composed of two sets X 
and U, 
ETS Tz <2. , n] 
U= [li,j)s i © X and je X). 


Pad 
Wt 








X is the set of nodes or vertices and the pair (i,j) e U is called an 


arc, where i is the initial node and j is the terminal node. Two nodes 
i and j are adjacent if they are distinct and there exists an arc u=(i,j) 
or v=(j,i). A path is a sequence (u;,u2, --+) of arcs of a network (X,U) 


such that the terminal node of each arc coincides with the initial node 


‘ ae 


of the succeeding arc. A finite path in which the initial node coincides 
with the terminal node is called a circuit. The length of a path is the 
sum of the lengths of the arcs forming the path. 

The matrix M associated with the network is a vertex-vertex matrix. 

_M= [mig]. | 

where the element mj; is the length of the arc from node i to adjacent 
node j, if no such arc exists then mj is equal to infinity. The dis- 
tance m5 can be less than, equal to, or greater than zero, but the 
length of any circuit must be non-negative. In general, mij does not 
necessarily equal M5 j and mij does not have to satisfy the triangular 


aS + 2 
ms J ms k Mk j . # 


The shortest distance matrix M* has elements mi j equal to the 
distance of the shortest path from i to j. Thus mj is the distance 
of the path from i to 7, whose sum of arcs is at a minimum. 

Following Murchland [15], the symbol ":=" shall denote "is 
replaced by." 

Following Bellman and Kalaba [1], 

min, (xX, Xo. tts Xp) = the nth smallest value of 

the quantities x,. = ; 


Symbols. frequently used are defined below. 


ND the number of nodes in the network. 


INF = infinity, and is taken to be greater than the maximum 


distance of any possible path. 
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All matrices in'the sequel represent vertex-vertex distance 


matrices. 


3. The Shortest Path by Matrix Manipulation 

The first method of determining the shortest path in a network, 
that of matrix manipulation [3], involves two operations. These opera- 
tions play the roles of elementary addition and multiplication 
respectively. 

Define the sum. of two matrices 

C=At+B. 

as 
b 


= min( 


ph (1) 


where aij and bij are corresponding elements. 


Sj Sigs 
Furthermore, define the product of two matrices 
C = AB 
as 


Cij = ASL + byj)- (2) 


M contains elements that are the arc distances between adjacent 
nodes. Following equation (2), mM? contains elements that are equal to 
the shortest distance, from i to j, in a path containing exactly two 
arcs. In general mk contains elements equal to the shortest distance, 
from i to j, in a path containing exactly k arcs. 

Thus, the shortest distance matrix, 

My =a Me Me 4 ee tyne), (3) 
contains elements equal to the minimum distance, from i to j, for all 
paths containing from one to n-1 arcs. 

As an example, consider the network in figure 1 and its associated 


initial distance matrix M, given below. 

















Flsiz 
m= | 1 {INF 4 
{inFl 1 LINE 





Using equation (2) we obtain for elements m2(1,1) and m@(1,2), 
m2(1,1) = min(7 + 7,5 + 1,3 + INF) = 6 
m2(1,2) = min(7 + 5,5 + INF,3 + 1) = 4. 
After computing the rest of the elements in similar fashion, we obtain 


the following matrix. 


9 6| 4 19 
Mo = 18) 5 | 4 
2 | INFI 5 











From equation (3) and using addition as defined in equation (1), we 


obtain the following shortest distance matrix. 


Mw = \6 
1 


2 








4 | 3 
a | 4 
Tis 








In genera] mn-l 45 a matrix using exactly n-1 arcs and visiting 
all n nodes. M" would be needed if a shortest circuit from a node back 
to itself was desired. In general 

Mo = M+ M" 
where M° shows the length of the shortest circuit from a node back to 


itself. For the above example 





5; WH: & 
Me | 6 | 5 | 9 
: oy & S 
and 
5; 413 
me =] | 5 | 4 
abt 15 








is the shortest circuit matrix. 
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Hetwork for Matrix Manipulation Example 





Fiqure 2 


Network for Cascade Algorithm Example 


1] 





4. The Cascade Algorithm 


Whereas the previous method needed n-1 matrix operations, i.e. ad 
for j = 2, + + +, n-] and one summation, the Cascade algorithm [7] needs 
only two matrix squarings. (This method is called the Cascade Algorithm 
by [7] and the Revised Matrix Algorithm by [12].) 

In this method newly calculated elements immediately replace 
existing elements. Also, now the order of calculation becomes important. 
In the first squaring, the so-called forward Cascade process, the elements 
are calculated in the order mj}, mj2, °° +s Mn3 M217, ° + + Mn; + + «3 
Ml» * + * Mn- Then in the second squaring, the so-called backward 
Cascade process, the elements are calculated jin the order myn; Mh .n-1> 

* Malt Melon» ° * * MPs? * <2 Mp ee ~ Mts 


Let MP = [mJ where 


: (} initially 
per | F at the end of the forward process 
_ at the end of the backward process 


Thus M! = M, which is the initial arc matrix previously defined 
and mM” is still the shortest distance matrix. 


MF 4s obtained from ml by the following forward process. 


F . q r 
15 + : 
ms mi (m, m4) (4) 
where 
ANT R23 
eee a9 
(5) 
and 
[ k= 4 
SUE k< i. 
a 
*e _M is then obtained from mF by the following backward process. 
Pies io r 
m, jrsimin (m’ +m.) (6) 


1 ‘ 


fj 
12 HH 





where _[f k <j 
; ata see 
and ' (7) 
IF oe | 
0 -(f k> i. 


To illustrate the above procedure, consider the network in figure 2 


and its associated initial distance matrix given below. 


014] 1 SINE 
ee PINEP OG] 1-f 1 
INF} 1 | 0 | 4 


INFIINFIINF| 0 


Beginning the forward process using equations (4) and (5), 


F as en 
m4 + mo =OQ0+4e24 


] ] 
Myo + Myyn = 44+ 0 = 4 
; my oi=min 12 id = Pic mM: 
m3 + M35 Sr++ P22 
] i L. 


Po Since 2 is less than M9» replacement would take place. 


mg:omin (0 + 1,2 + 1,1 + O,INF + INF) = } 
. Since m3 = m3 no replacement need be done, although when using the 
computer, replacement would automatically take place. 
mygi=min (0 + INF,2 + 1,1 + 4,INF +0) = 3.< my 


Therefore mat 4 would in fact take place. After computing the remain- 


ing elements and replacing when necessary, we obtain the following matrix. 


Op oy t Fa 
dF a PINE) Ged ft 
! INF] 1] 0 | 2 


INFUINFJ INF | O 
After performing the backward process using equations (6) and (7), we 


obtain the shortest distance matrix M* as shown. 


. epee 12 
ues MEL OTT Va 

_ INF} 1} 0 [2 
INF LINE! INF 1 0 





Thus we see that this is exactly the same as matrix squaring with 


‘ the exception of replacing the elements once shorter paths are found. 
This replacing cuts down the number of operations performed tremendously, 
relative to the matrix manipulation method. 
Narahari Pandit [16] first claimed that two forward processes were 2 
enough. However, [7] displayed a counterexample and proved one forward 
and one backward process was sufficient. Hu [12] later proved that Shves 


forward processes are sufficient. 


5. An Inductive Approach 
Dantzig [6] takes somewhat of a different approach to solve the 
shortest path problem. His inductive approach is described below. 
Assume for nodes 1, 2, - + +, k-1 that optimal distances ms are 
given. Optimal distances mi are desired for nodes 1, 2, °° *, k. 


For h=1,°°*,k-lTandje=l,-°+sk- 1. 


* 
¥ Meh = min (m5 + mh) _ (8) - 

‘ = j . oO | 

"hk min (m5, ; hs) _ (9) 

m = min (0, Mme +f Ms ) (10) 
kk ny kj © “5k | 

For i dye « og RS Vand J =), + ey kT 

* P * * : 

my gs=min! (mg st x + m5) cc (11) 


Consider the network in figure 3 and its associated initial distance 


matrix given below. 


0} 1 |5 {INF 

ge > JINEL OD PF] 
INF| 5 |0]} INF 

: INFIINE 111 0 


Assume we have found optimal distances ms for nodes 1, 2, 3 thus 


giving Mo. 


O;1{ 4 
1 | Q| 3 
5 + 0 
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Figure 3 


Network for Inductive Example 





Figure 4 


Network for Algorithm 97 Example 





: * 
To obtain mij for nodes 1, 2, °° *; k,; we begin by using 


equation (8). 


! 
Ma, * ™h 
* : 0 
Mg, = min Mao + Mp 
j 0 
Mg3 * Mp 
For h = 1,2,3 
* 
mg, = min (INF + O,INF + INF,1 + INF) = INF 
* : 
Myo = min (INF + 1,INF + 0,1 + 5) = 6 
43 = min (INE + 4),INF + 3,1 + 0) = 


Concluding the calculations and replacements using equations (9), (10) 


* 
and (11) we obtain M . 


Cart Teh 2 
oh MEE PEL 

INF{5 | 0| 6 

I INF!6 1110 


6: ‘Algorithm 97 


Unnoticed for some time was the following algorithm by Floyd [8]. 


The two «'d statements were not in the original program but have been 


added to reduce operations. The program is written in FORTRAN IV. 


DO 1 42.1, ND 

DO | d= 15 ND 
IF(M(J,1).EQ. IN \F) GO TO 1 
*IF(I.EQ.0) GO TO 1 

DO 1 K-#i, ND 
IF(M(I,K).EQ.INF) GO TO 1 
*IF(I.EQ.K) GO TO 1 
MS = M(J,I) + M(I,K) 
IF(MS.GE.M(J,K)) GO TO 1 
M(J.K) = MS 
] CONTINUE 


Flow diagram for the above program is contained in appendix I. 


tion as the Gas cede algorithm, 
My? =m Mey ss 4 + mp) 
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This is similar to the previous matrix squaring methods. It is not 


actually a matrix squaring but resembles it. It performs the same opera- 


be‘worked to illustrate the procedure. 
non-infinity, non-main-diagonal elements. 
sum is less than the third element. 


The search is as follows. 


searched in order, k= 1, °° + 





but: calculates Mey in a different order. Thus, by this order needs only 


a single "squaring." Upon "squaring" Mik is now Max: 


After a brief informal discussion of the algorithm an example will 


The matrix M is searched, in an order to be described shortly, for 


Two of these elements are 


added and compared with a third element. Replacement is then made if the 


The columns of M are searched in order, 
Each column i is searched in order, j = 1, °°: > 
When a non-infinity element (j, i) is encountered in column i, row i is 


For all non-infinity elements in 


If ms < Migs then M5, 2=Ms. Thus we have formed a shorter path from j to 


.k by combining two paths j to i and i to k. 


distance matrix. 


——O~ MW 


INF INF] INF | INF| INF 
INF{INF)INFJINF) 2 
INFtINFIINFtINF) INF 


column 1 has no non-infinity elements. 


first comvarison, 


As an example consider figure 4 and the following associated initial 


INF 
INF 
2 
] 
INF 
0 
INF 
INF 


Applying the algorithm to M we see that except for i 


Incrementing i to 2, the first 


non-infinity column element encountered is mj; = mj2. Now row 2 is 


searched and the first non-infinity element is ms, = mg3. 


Since 2 = ms = m3 = Z, 


Thus for the 


Re ITT 





no replacement is needed. The two remaining non-infinity elements in 


row 2 are the next elements utilized. 
Se NS = Mo + M4 
b4.2 23 


3= ms =m, = 3 ~ 
mgt? ms 
3. ise = Mio “ Mor 
Et bee 


nm 
ul 


ms < M5 = _ INF 
My 5i=2 


For i = 3 the following two replacements would take place in order 


aiven. 
4. mg=4 
= 5. Mog?=3 + 
After the remainino columns, i = 4, ° ‘ *, 8, have been searched 


and required replacements made, the shortest distance matrix is obtained 


and is given below. 


wal OGptte¢S5 |e) 416 42 
INF} OTT TALI STS 16 
INFIINFiO}1!16121/4 45 i 
INFIINFIINF}O |} 5/1413 /4 
INFIINF{INF} 1} 0} 214 |5 
INF/INFIINF} 5 | 4] 0}]2 |3 
INFIINFIINF!} 3 | 2} 41]0 41 
INF }INFIINF LINFIINFIINFIINE | 0 


Floyd [8] gave no proof of his algorithm but celia for proof ona 
theorem on boolean matrices by Warshal] [19]. Hu [11] gives a somewhat 
simpler proof of Floyd's algorithm. Murchland [15] developed a similar , 
” alaorithm, unaware of Floyd's accomplishment, and presents a somewhat 


lenathly proof in his paper. 
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hd 





If we are interested in not only the distance of the shortest path 


but the route as well, something more than just M* is needed. One way 
to obtain the path is to define a new matrix, called the routing matrix, 
MR =[mr sy] : 

where MPs contains the number of the first node, after j, on the path 
from j to k. Initially, | 

mr sy =k, FOP ul« KS Lew’ * petits 
During the calculation; whenever 

ln i Se 
the following replacement is also made, 

mrs, iamr as 
In the previous example the first replacement is done in comparison 3. 


Thus 
@ 5 


Ba 
would be followed by 
Mrs t=Mryo = —— 
At the conclusion of the above example the routing matrix would be 


as shown as follows. 


MR=]17 [2 [3 [4 2 [3 |4 | 4 
i2 jo 4 [e481 a re 
Yy2ei3 (4 /@ 1 @& 7a | 4 
Ts| 2*{ 34/4 12 | £27 Fs 
1% 1 23(/ St; 4 |& tL 1? 74 
Tegey ie tee Te RS 
A) 2") Si StS 1S 17 Pe 
Per ab ST ae Poe oe hee ae 


The elements with the stars are elements as they were initialized. 
No replacement took place since no finite path exists, as can be seen 
in M, In the computer program these starred elements would be set to 


infinity to facilitate reading MR. 
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To illustrate the procedure of extracting a path from MR, the path 


from node 1 to node 8 would be 


mr = 4 
Mag = 8 

or 
i to 4 tor &. 


It may be observed from the example that more than one shortest 
path exists from 1-to 8. Algorithm 97 will only give the first best path. 
The modification required to extract all alternate paths shall be treated 
in the following section. 

The above procedure of using the routina matrix to determine the 
paths is applicable, not only to Algorithm 97 but to all four matrix 


methods. 


7. Comparison of Methods 

Already some idea exists of the relative efficiency of the four 
methods discussed.-- Simply in terms of number of matrix squarinas, the 
matrix manipulation method requires n-1, the Cascade algorithm two, and 
Algorithm 97 one. To more critically compare all four methods let's con- 
sider the basic operations involved. A matrix multiplication between two 
matrices of order n requires three basic operations. These — 
are indexing, addition, and comparison. Below is a listing of the four 


methods and next to each, the number of each of the basic operations 


required, 
Matrix Manipulation - n? log (n-1) 
Cascade Algorithm - 2n° 
Induction Method - n(n-1) (n-2) 
Algorithm 97 - n(n-1) (n-2) 
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8. All Alternate Best Paths 


In practical problems it may arise that all alternate best paths are 
wanted. A decision maker may want to evaluate all best paths and select 
from them by means of some other criterion. To extract all best paths 
Algorithm 97 must be modified slightly. 

For the shortest path problem we needed two ND x ND matrices, M and 
MR. For this problem M will remain the same. Recall that MMs contains 
the number of the first node, after j, on the path from j to k. Now 
since we want to keep all best paths. we must go to a three-dimensional 
matrix. The third dimension on MR is not governed by the number of best 
paths from j to k but by the number of first nodes on the path from j to 
k. A number of paths may share the same first node and later branch off. 
Consider the paths u = (1-2-4-7-8) and v = (1-2-4-6-8). Node 2 need only 
be stored once, since paths u and v coincide from node 1 to node 4. Indi- 
cation of a branch would occur at node 4. Hare two first nodes would 
appear, node 6 and node 7. 

A priori it is-not- known: how many best paths there may be from j to 
k. Thus it also is not known how many first nodes are needed. This is 
of no real concern if calculations are qoing to be done with pencil and 
paper. For the computer, however, a specific number is needed. Inspec- 
tion of the given network can give an indication of the maximum number of 
first nodes expected. Let NA denote this maximum number. Surely 
| NA < ND. 

Thus MR would be dimensioned ND x ND x NA and mM sg would contain the 
‘mer of the ith first node, after j, on a best path from j to k. In 
the computer program, to facilitate printing out the paths and readina 
the getriees: mr skin is set to infinity for h = 2, * © °, HA. Also, upon 


completion of the algorithm, for all elements Mek equal to infinity, the 


2] 


corresponding elements mrp, are set to infinity. Thus the routing matrix 


will also indicate all non-existing paths. 

For the computer program in appendix III, a value was read in for NA. 
If during execution of the program the number of first nodes becomes 
equal to NA, an aviar message will be printed and execution will cease. 
This will prevent the computer from attempting to address either outside 
the matrix or wrong elements in the matrix. Since computer time may be 
costly the program could be altered so as to ignore all additional first 
nodes once the total number exceeded NA. An error message would still be 
desired to indicate that alternate paths may have been thrown away. 

Now the modification of Alaorithm 97 will be explained to enable us 
to extract all best- paths. Departure from the main algorithm occurs at 
the point where path (j,k) is compared with the sum of the two paths 
(j,i) and (i,k). If the: sum jis: greater, that path, as before, is ignored. 
Now let us consider: separately the cases where the sum is less than and 
where it is equal to ie 

For the case "less than," a path has been found that is shorter than 
the existing path. This is sie ware as for the shortest path. Now, 


however, 


mr. pel,” * *, NA 


ike Sip 
In the previous section for the shortest path, p was equal to one. Now 
there may already- exist more than one first node at mM sips 
all must be transferred to Mr skp* (It is true that there may not be NA 


first nodes at mrs The transfer could be handled a number of ways on 


ip’ 
the computer. Qne way would be to determine which of the NA elements are 
first nodes and transfer only them. The other way is as above, transfer 
all NA elements. Whichever is more efficient would depend on the par- 


ticular matrix.) 
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Therefore they 





sae Sibvoenirisains 


eA, tah athens 





For the case "equal to," a path has been found that is the same 


distance as the existing path. This we want to keep. Now a comparison 


must be made between: the first nodes of M skp and mr Only first 


Jiq’ 
nodes of mM sig that are distinct from the first nodes Ss Mr skp will be 
transferred * mM skp Redundancy of first nodes would only take up 

| storage space and no paths will be thrown away by the discarding. 

As an example consider the network in eure 4, M and M* remain 


unchanged. The routing matrix would now be as given below. 


MR = || 1 2 3 fe A ]253 |) 2 3 12{4]| 4-]2)3 | 4)2)3 
INF z 3 A 19151 & 3 {4/51 4 (3/5 1 413)5 
INF INF 3 4 4 \6 6 14 4 16 416 
INF INF INF 4 7 |6 6 7 16 816} 7 
INF INF INF 4. Ss, a 7 \4 4\7 
INF iN INF 7 7 6 q 817 
INF INF INF a 5 5 7 8 
INF INF INFI. INF INF INF INF 8 





The columns of MR are: divided by double lines. Between the double 
lines are the elements of the third dimension, the first, second and third 
first nodes of the alternate paths from j to k. It can be shown that 
there are: 25 alternate paths from node 1 to node 8. 

Appendix II shows the flow diagram. for the algorithm. The full 


computer program is shown in appendix III. 


9. Nth Best Path 

A cdpepal zation of the shortest path problem is the nth best path 
problem. It may arise that-the shortest path(s) may not be the "best" 
path when evaluated by some other criterion. A decision maker may be 
willing to deviate from: the shortest: path(s)- by a certain amount to avoid 
undesirable aspects of the shortest path(s). Thus the nth best path 
becomes important. “ 


First a brief review: of the known methods is in order. The simplest 


method is given by Pollack [17]. .The network must first be solved for 
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the shortest path, say u. Actually all shortest paths, Up? must be 


found. Once given Up? from a given node j to a given node k, composed 
of ry arcs, the length of each arc in Uy is set, in turn, to infinity. 
The shortest path problem is solved (at most) ry times for each Uys once 
for each arc set to infinity. (An arc appearing in more than one Uy 
need be set to infinity only once.) Thus the maximum number of shortest 


path problems to be solved is equal to the sum of the ry Ss The short- 


est of these new best: paths is the second best path, say v. There may 


- be more than one second best path. Thus Yq will denote the qth second 


best path and will have o5 arcs. 


’ For the third best path, the length of each arc of u. and u. is set, 


Pp 
in turn, to infinity. Extending this to the nth best path, the length 


of each arc of all first best, second best, * ° *, n-1 best paths must 


-in turn be set to infinity. As before, for each arc set to infinity a 


shortest path problem must be solved. By now it must be apparent that 
the number of shortest: path problems to be solved may quickly become 
astronomical ly large.--Thus this method appears: to be good only when r 
and s are small and when p and: q:are very small,- hopefully one. 

The method of Bock, Kanter, and Haynes [5] differs from Pollack's 
in that the shortest path need not be known a priori. Using stems and 
trees, this method simply is a systematic listing of all paths from a 


given j to a given k.- Again, however, this method is limited to smal] 


networks. 


As in the two previous, the method of Hoffman and Pavley [10] 
determines the nth best path from a given j to a given k. However, now 
not only the kia ledge: of the shortest path from j tok is required, but 
the knowledge of all shortest paths from j to all other nodes in the 


network. Then by deviations from all of the shortest paths the nth best 
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path is determined. Thus we have a problem similar to the one in 


Pollack's method. Pollack's method, however, being a shortest path 


problem generates paths with no loops. Such is not the case with this 


method. Paths with loops may well emerge as nth best paths. How much 
of a problem this is depends on the particular results desired. 
The limitation to small networks js escaped by Bellman and Kalaba's 


[1] method. Here the nth best path is determined from all nodes j of 





the network: to a given node k. This can be considered an advantage or 
disadvantage depending on what paths are wanted. This method uses the 
functional equation technique of dynamic programming. Restrictions are 
imposed in the method to insure loopless paths. 

Floyd's algorithm provides a basis for the development of the 
following matrix method to determine the nth best path. The method of 
searching the matrix remains the same but now additional calculations are 
performed. Upon completion of this method the first, second, * ° *, nth 
best paths are given from all nodes j, to all nodes k. The only restric- 
tions placed on the size of the network, or how large n-can be, is the 
amount of computer storage available. The description of the method wil] 
be given below with the flow diagram given in appendix IV and the full 
computer program in appendix V. 

For the method to: determine all alternate best paths the matrix M | 
was unchanged, but: MR was increased to a three-dimensional matrix. To 
accommodate the nth: best: path both matrices must be three-dimensional. 

M and MR will be dimensioned ND x ND x NK. NK is equal to the pre- 
determined value of n. Now, Mako will contain the distance of the pth 
best path from j to k and MP skp will contain the first node of the pth 


best path from j to k. 





To determine the nth best path from j to k, the following operation 


is performed. 


Mikn “7 mintM sen? min, (mosp = Mig) 
a Ty * * eRe : 
a2 ty 9? swe 


Essentially this operation forms all combinations of paths from j to k 
through i. Since not all of the NK p's or q's may be non-infinity, there 


2 


will be a maximum of (NK)° of these paths. The path through i with the 


nth best distance is compared with the existing path, Mekn The minimum 
* * 

of these two paths is Ms kn The Mkn replaced by Mtn 18 not simply 

discarded but becomes the Ms ne best distance. (It can be observed 

that when n, p and q-are all one, the above operation reduces to the 

same operation originally described in the Cascade Algorithm and 

Algorithm 97.) 

To illustrate the procedure consider the network in figure 5 and 
its associated initial distance matrix given below. 
M =) 0 2 4 9 | INF 
INF | 0. | 6 | INF] 7 
INF | INF} O 4 | INE * 
INF | INF | INF} O 8 
INF |} INF 9 ! 13 0 
The first non-infinity column element encountered is Mein = M101: 
Searching row 2, the first non-infinity element jis Mery = Mo31°- Follow- 
ing is the first set of comparisons. 
1.1 ms = M91 + M3) = + & = 2 > gy 7 4 
M134 :# ms 
2 We= Myo, * M534 =. 9 < M139 = INF 


M1 39 2= 9 
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Network for Nth Best Path Example 


a 











Continuing across row 2, 


ZL. te > M194 “ Moe] = 3+72+10< M5] = INF 
M57 = 1G. ; 
Following are the steps for column 3. 


3.1 ms = M31 + Ma] = 4+42=8< Ma] = 9 7 


Wt 


mh de + 
| Taal * 
3.2 HS = M139 + Ma, = 9+4= 13> Mi 40 = 9 
M40 :# ms 
3.3 ms = Myo + M3q) = 13 < Ma = INF 





M43 °* 13 
4,1 ms = Mp3] + M3a] = 6+4+ 10< Moay = INF 
Moqy °* 10 
5.1 ms = Mp 37°+ Mga) = 9+ 4 13 = Men] 
Mey 27% ms 
At this point the matrix would be as given below. 
M= |} 0 3: 4 191} 1} 8] 9} 13|]10 | 
INF 0 6 10 7 
INF INF 0 4 INF 
INF INF INF 0 8 
INF INF 9 13 0 





Upon completion of the algorithm, M* would appear as below. 


M* =|1 0 3 4) 9119) 8| 9413{|10 116|17 
INF 6116}271)10|20131|| 7118 
INF 0/21 4|25 12] 
INF 7 0/21 8 
INFL |} I 9 ie 0121 


As before, the columns of the matrices are divided by double lines. 
Between the double lines: are the three elements of the third dimension, | 
the first, second and third best distances from j to.k. 


Proof of the nth best .path algorithm relies on the proof of 


: 
i 


Algorithm 97. For if Algorithm 97 is valid, it must be true that all 
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possible combinations of paths are enumerated and only the shortest 


retained. If Algorithm 97 fails to consider even one path from j to k, 
then no claim can be made of its validity. For there can be no guarantee 
that this one neglected path is not the shortest. The proof of Algorithm 


97 has been sufficiently established by Warshal] [19], Murchland [15] and 


— Hu [11]. Thus,-all possible combinations of paths from j to k are con- 


sidered and only the shortest kept. For the nth best path algorithm, 


the remaining: paths are not discarded. All paths are essentially placed 


_ in order and stored. 


The procedure employed in constructing MR to enable the path to be 
enumerated is the same as before. First node replacement occurs immedi- 
ately after the distance replacement. Considering ‘i preceding example, 
replacements in MR would occur as follows. 

2 Mry35 °= Mryo_ = 2 
‘2.1 Mra, °= M407 = 2 
= 4 


Ba mig. : my 4] 


miiap yg Pe ie 
3.3 Mry a3 $= MP1 99 = 2 
At this point the first row of MR would appear as the following. 
lay a det t wsra1 wstateauert | : 


Following is MR as it exists upon completion of the algorithm. 


MR = Ih or 1 uzi2}2 isi 4t 221314 
2 3/513 13)5/ 315] 3 
3/4 A\ 4 4 
5 A| 5 5 
a 4 514 


Extraction of the path from MR is somewhat different than in the 
previous examples. The first node of the nth best path from j to k will 
not necessarily be found in the nth, third dimensional element. Con- 


sider the nth best path u= (j, °°, Pp» ° “s " *) kj). The first node 


29 


of the path u would be mr and found in the nth element. Succeeding 


. jkn 
first nodes may be found in the nth element until say the pth node is 
reached. Now the first node of the path from p to k may be in the n-Ist 
element, indicating that this is now the n-1lst best path from p to k. 
Continuing, the path from q to k may be the first best path. Thus, the — 
distance matrices (including original) must be used in conjunction with 
the routing matrices as the path is traced from j to k. Upon reaching 
the first node of u, Mr sin? the distance of the arc from j to Ms kn must 
be subtracted from the total distance of u. The elements of TMiaedt are 
now searched for this lesser distance. The third-dimensional element 
will iuateate what best path we shall continue on, 

Using matrices M, M, and MR let us determine the second best oake 
from node 1 to node we From mry55 we obtain the first node as 3. Sub- 
tracting the original M134 from the total distance of the path M150 = 16, 
we obtain M55 = 12. Searching Map. We find i = 1, meaning from now on 
we are on the first best path so no more subtracting is necessary.’ Con- 
tinuing, Mr351 = A and lastly Mas, = 5, Thus, the second best path from 


- 


1 to 5 is u= (1, 3, 4, 5) and has a distance of 16 units. 


10. The Traveling Salesman Problem 

Algorithm 97 determines the shortest path from j to k. The first 
extension discussed allowed us to expose all possible minimum paths. In 
both cases ms was initialized to zero. If, instead, Ws was set to 
infinity, we could have extracted the minimum circuits. Extending this, 
way m5 was set to infinity in the nth best path algorithm, we would fe 
obtained the nth best circuit from j to j, for all j. If the nth best 
circuit: contains-all the nodes of the network and the anty node appearing 


more than once is j,-which appears exactly twice, first and last, that 
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Se 


circuit is defined as a tour and is a contender for the optimal solution 


to the traveling Salesman problem, Thus , the traveling salesman problem 
modeled-as a nth best path problem js next to be investigated. 

The nth best path formulation in the previous section is not a 
minimum path formulation, as is Pollack's. Thus, paths with loops can 
be expected. This could be remedied by the insertion into the computer 
program of a subroutine called, say, NOLOPS. Each time a new path was 
formed by adding two paths, NOLOPS would construct the path and discard 
it if it obtained a loop. 

Thus we are assured of loopless paths. Recall that alternate nth 
best paths were discarded. Only the first nth best path was retained. 
There is no way of knowing if that discarded alternate nth best path was 
optimal. A way of fixing this would be to compare alternate paths as 
they arise and retain: the one with the most nodes. But this is not 
really correct. For that alternate nth best path thrown away may have 
combined with some other sath and formed an optimal circuit. Thus, all 
alternate nth best paths must be retained. Now, all alternate nth best 
combinations of paths must be considered. This could be accomplished by 
adding a forth-dimension to the matrices M and MR.. Each would be dimen- 
sioned: ND x ND x NK x NA and Ms knp would be the pth alternate nth best 
path from j to k. Upon completion of the algorithm the n best circuits 
from j to j would be searched and the first tour encountered would be 
an optimal solution. 

These modifications would seem to overcome the difficulties of the 
nth — formulation, The optimal solution would be an exact solution to 
the traveling salesman: problem. It is anticipated, however, that to handle 
large networks in this way a computer would have to satisfy the following 


conditions. 














Handle four dimensions 
Have a storage capacity larger than even the present 
day advanced computers 


Be faster than the present day computers 
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APPENDIX I 


Flow Diagram for the Shortest Path Algorithm 
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APPENDIX IT 


Flow Diacram for the All Alternate Path f° gorithm 


CION 0 oa 


ay a a 


FALSE 


’ 


MS=M(J,1)+M(1,K) 


> 


MS-M(J,K 


< 
M(J,K)=MS 
15 
L=1,NA 


5 MR(J,K,L)=MR(J,1,L) 

















JOVYOLS JTGWIIVAY Sd3d9x4 
SA3GON -LSYI4 JO YaaWnN 


) yin 
4NT=(N°X°C) UW } 


6 $a Caer ) 


Wi) 
an Bl 















sa ac ah et ie me he aE ai ik ch a aca ete : 
tt tt el Se tl ct ec at 


p OUT 


APPENDIX ITI 
COMPLITER PROGRAM FOR THE ALL ALTERNATE 


MATRICES TC RE 


LTERNATE ROUTING 





LATE PRINT OUT 


1,NN),1 


38 


GO TN 29 
GO TO 20 


MATRIX 


tei 


CaalseCn a 
, ASMIUOS ee 


DEROREXU XC SEL 


NISTANCE MATRIX 


ee Be 


i “OFAC eEZ 
rite SOwuw Liste me ZONA ZS 

: Me YOO mix ZO. 
MS eu eh Zea ZOW 
a2 BO rent elie 


NITY ELEMENTS OF M(I,J) 


a es 


AL MATRIX TO SHORTEST DISTANCE MATRIX AND 


e(atisedbédt 


NUTING 
FQ.INF) 


i) MWD CWC OWS 7 Nw AAU et NZ 


a Soe = 





Beene e ec 


OSS SSS (STS SSS SS SiS) 


LYS RORY ROR NI fat SHO Zee HOWR SE It ile ell eomeme Ss OT wl om eu 
NO BMNAOMZ Line eM of Oe MOm~SPOZe es etry 
ae wt et ere oF LY 


“OO eS Ho eNIMZIMD 
Tl mIUQ in 


e eH ff eM PWWROZWOH AAA & PZOOZZAZEFMONNE HONS REE eA ORAS AT Se 
22 Ke IW Ser TOSS SIAM UOr 
pet HLL | HO Hem MO 


i IO TO 
BeNOR Dt WYETH COC CK e™ WRK KYOYOWTCOMUTUUMALK-OYODOLOWUO 
ZZ VORORY eR Dis ZOOZALCA HZOOZVOZBDZEOCAAA OOP LOD HH Ores ae TOSZOC NOM 
LUNA IS dA OMND~ aw 


were eee [fe De 





ae, 
re = 
& Lu 
oO OQ 
a << 
Cc [24 
Pa eu) e 
pan & _ é 
eS a vg) ; 
oO ~ SH [ 
Qo. = ~~ SO j 
Ww c ~ ew ‘ 
tu \ aw o we tu } 
eg ue «x «KO 
(a4 — Oe 
v2) eS 4) cr ew 
uw oO — — bk 
¢ be ~ <M, t 
Cc ea) <q Ere | 
=z 4 i G: ~ ~~ OC 
< >< = - ~WS OG 
ke = — * ~ Ww NO eZ Ze 
Vv) »~ [ag al ~ ey e ow tw \ 
a ol Daal oad be lead ” - < <a oo IKE 
ci + ee << ~ - = - = > be ee DU 
UL. << ms om = ~ & a N ec —AY Oat 
”n e 4 ze ian fo! Pas ond z Ww OL th OL 
tis Uh it mL Ww = Ls " i al Rot Ww 
Oo O <— we G Oo * N x pas =) < Bw 
met he Ze ee Zz a < o ~ ain Se RIL ~ a 
fa dian yaa = Hf f u hoo ~ - ina) MEaxO ~ aN 
re Ww — im Oo ke 2? f= 3D a Dm es O 4 wUIZZ Oo OWN 
-~ TON w onl ga YA - = NO ~ ~ = NO ~ m7 a © > eal eel ® oS a mt 
- a aa. a Aa Zz -~ & ma Mw Zz Mw << k- tO eR LU ree 
. = — ~~ - a 27" oe + = = & Os « & FODrRO “SAN StH 
aiid FOZ >. if —e a - e RN - = ee ) EN 7 ke © STO Sees 
- a5 re  <CO eee Zz re a Zz a z~ = Zu xrgtD aO ann) 
3 oS = a et aN _ ww tr - On ~~ te ore fs © oO _ . =z NOTA ® 
— ~ e lad ye ee OS ee Go ae a = Ge) S&S ww ~ ere a a oP a 
ce Foz CaO Ase — WwTZOOr!Z &F “QO “O—- = Tt »O~ -Q— ~- © - & & of O * LOC LN 
= Z2emw 22 260@ 42 OC 22a XY SZ ew ewes ee He ee me oz - ee LA MAH & 
If ww TO ~ ns = Ye ew ey Oo wee Ow ee DZD ow © Owe mm OO Aaa HAO & & & ord 
~ OW<l set et © NOetaem LT ertae eorte CO Fert eee 8 YS Wi CCITT oe KMKnL 
MNZWUO eS It tt ml be Ur OUI et HNL HNL I Tm NLU) J] NLU DOW Om Wet tra Or Otte 
an et 3 Me TOUDAYeg & ed i neh ar Nada ae a ee ZeifjDe eS eee 
Xe Zzadi a a ee Sa ae Pl ia © 42qnW IGA Zee RRR Ree eee 


rO eR EZ NN Fae Tae SMOe bea ae aa an a feo oq weeded 
Zi RR ZR ANNONE YR WWNN EHR OZZNZ RE ZSIZNZREOZMANZ RZZAN 2 eer Zee ZH TEST TTT STTSAS 
ee KS I fT ee be oe a PP Pa ee on Pe ee ames wade ae Bede a an) 
SYDYODZUYTOOTOVOLOWWOOW YOR LYTYKOULOYOFYOYVOYOUCOYOLYOAOLANOTIOO0OVO0C0ANoN2Z 
WOSZQUZ WOO YOO Fc aA Na ta a i ac la lat 
Lead Ot LL _ 
MOow a NO mat + Vale ae) - DR OW FB NOPHNOEDOOANM 
aA O N ALL Na N NON N AN Me OH MAMMMOMMO Sry 


oO vo QO i : Oo QO 











eaten a —_— a ar ne neu nee ian tlt nl ipa a a in i onl at ent LD a eT a 


TO NODE', 


ce 
on) 
tamed 
td 
uJ 
a 
Qo 
= 
= 
Oo 
fo 
fam) 
za & Wi =m 
~ bk WO ow» 
=O = WM” 
Ys Cc nN < x 
coe ae co a Enel hk 
ZO + ~ + wn < 
WILY -~ eae a me oO to OL oO 
zu LL = —_ - <a Ge <= 
x 2 Zz C = . + = WE 
Ot _ ke a I a Cc co i =D 
Or FF ° = = wo © —- ‘= oo ie 
ae ~ Ww (=) oO <I he mi con <I bm» b> al] 
z2zeOS G@ & uVoo _ =z <a o ~ Lod = < me 
rome ~ Ga) rg e co = —O = Za et 6 
We eee ac o“a (ou -_ Co wmf ° je = o me = 
ee -OOZ “+ cic oo m & - ao [e) wa 8 iw > 
AOWZ = eee ~e@ ee Zz aD ui e z= a Me MO 
<2~wWwO {a ed ULL ~ Go w —OO e ~ OC W ~ = eet 
eZZ0090 CO AtoZ & a ad IZ a OW Ze = 2 a 2H Ze + ee © 
Luu eZ 2 2wWN Zw o~atZ me Ze eX wm ert Ut me ee OX ~ te QO) 
ZZZO « o z= azz Onn & xa aoe SELL eOrt } Z2ZeawdW wie SW + TON « 
SOOM OaNOe eF UW aay iY Yt OY LL ee Delf — LRH eu tl Le ee a We 
eRe Steel ft ie wenn Ws eH noe fF eDatit : Seite ew DOW Ne et eDoe II a-Zzto_~ ~~ 
DZVi~w ae heokent A etme ONIN mer OOO oh atid et Kae ee MZEa~ FPO Zee ers 


COZERRE — hee CO BHO pew AZ Edt SOK AO Le (eA SHE HM |] VAX | MOMS odie 
CT WYSSY eH oe WOM an e¢-KeESOt+ © WEP ODER Wee VA ZHS 1 oe oj HEAR oF Weer EMEED 
OTe A weNZOHZ AMADRenH ff ZO wes ew ZK HZ lew ROR ON Zee eZ KAM ORL OO 


DOSY LO~OU~ OH Ow HULL ZO HS NH UO NOLO fie On WU TOUOWN Tt WL TOOD OCOwe 
NVOZLACEOHSZOZVXOL2 EAH ZOZz BEM EZOHOZE MMOH ANOMGVZODPZEZTOMOZEVI~HIOOW Witaclu 
od 
AN ina] Baila) of amo is) aN ot Uy OP OO 


-7 ee ae ee eee 





APPENDIX IV 


Flow Diagram for the Nth Best Path Algorithm 
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ih: 
ey oat (> LL} Qa 
o eZ - = 
a8 Zz = (22) a 
peat <t bad uu 
oO w rt Qe foe) © >< 
roo MW wo x oa = 
<r unc i uw 
O CO LkO ~ >< Ww Cc ad 
QO —- Sew OQ ALUN = 2 < 
ae Ze ta) NZ om CG = 
nO Onn - Ze>O be a. 
uw uw <ttiL lw | ek 2 erie =) [Vy] ie) 
oO oO za pera.) ae Q uu W 
Oo =m om oOo~ ZeZze ~ of z 
RPRaZ2 aw 7D =F aad Co eed OO Ss & < 
Se) - Ore b= eel 4 = = (rae Lael 
Caz 2a eCc~y ae qeZz = Oo ” 
Oo Ft k&e ro =z Zeal —_ 
eee WN one ae & me Zz LL Qa 
Wut CD mle foo Pm eroa — Rw 
ke pa Ce of Lee = ‘ «I - 
ZZW = Fee wom ey FT Lu < _ vy 
em Te C$ Df ELAS ~ZreO ib ~ - ~ Wd 
Lee TS YO Tee Zz euiex eS = ag ~~ $e t+ fea] 2 (35) 
aaq Cc & aos elu - QO ee a av] A! 
a w > eae mel az Mm a mete a 
tui MOuws Ora Zul & - = . ae o Land cS. ey. 
Ode Fe WY ats W2ZXs = WL a Rl. Lonel x b— Land 
Mm MW OUZ LF Ww Zz et io Ly Ww "i u. Zz. x< 
noow w Zz = eZ NeD Mi = = S) 7 Zr oO Ou () Q 
Teeo BS HOt KOM COU NN ve = ~ wt} o bo oo [ae] 
-— Fi Reo =r NO WOZYre — hal a = t-~ = 
ZIIT¥Y Y¥eM=A eC AO PZ & HOU us oe ee -—~— At ee me << OQ -O© 
CO. Om YW FO eZ W e¥Y MUS ~ oH & »- “eo u as WAN LLN 
Iq wzZzoud <wu%C + EYZe et >= ee — ££ =—3 =. x Ze z 
roauw — ZPWMM SIZ ME ke = aD ~ -— _ reo sO HO 
nw aI fuuoes DT eZrne -O—b~ a a « diz Ce 
LUO UL BZN 0 OO et AOZES]SY a & ol y — Be to > a. Go 
Qmocw WwW abe O2mnine Dag e ie WL tee et Xt a = af i e wo wd 
oO Me LOK Ah me ee oO. euwust OO WYWLSZ Le OC oe PO~ LFEOO eZ UTC eH eD 
ZMUWWAO oe RUE MZ ZOWNZ ZH ZZ ZH ZDeSZere mee O ZZznn DOZZan Sm 
DONWZWHEAULXHD RH WRATH Ut ww Ol oe |} OM & Owe S or Tn bee a 
K~HLULODW WOROW OWZZ AR MMW RAINY et | Zens Ss eet ® DOs we ted eet Oo 
eCCZZOuSOd SKE OCHRY OR ee [LR ae~R Zar He SN OI De tH oo 
_ “i SIIOQD ee NAOMZIMHD AMY eAO © BHR NMMD ZY © wD ZI D “OGY «+O 
ZaeterlL LOPOUAaAsaeMZsze MM Wa — o er ees bar) oad Zui eZ “Uy tL 
KWON SHOm Ue SEMA ZOIO RPE AMIAO OQ wee we NG SDR RMR eS Sw er UOC 8m 6 
urord ~ -Yrat MCUEWZOZUID ata eet A eZO ET AR HOZ ZA ZR WW HE ONE ROIS 
SEERA IAORZ ee LIMITS SaaS St et ee bet et Oe Pe be wow ZK ed 


HSS NZNO en_neOe WX ROOK HK OV eS WWOOMR Yr eCOXCOUWOOUWYOWdIOOULOUL 
WZZNODWEY eH DES ZOVZOMKA KH KCOOPDOCZARCHAMOOZZAACAO BAQORMZOS> COC ee One 
Mt HOSeZtaAHOQnn~ oe <— as. <i GC) — oc ly Ze 
Zzaxr~ouwr = ~ [a a] us = Ort Na moO TO 
ezzzzz Zz =F =z “ <4 aay YN AO pire AO 


VQVOIUVOQOVOVOOCOYUOUH =) OOU WY oO oO STS) VO 


43 





ke eet pnts ler gt 


G 
- 
Qa 
[ao 
a 
Ww ne 
a ns 
(an) os 
Pam -~ + 
<= 
= mi 
GC [= ahaa! 
[a4 al * 4 
Lo at 
=< 
— =z 
beret ae wm LJ 
” Late) Ke = - ™~OU 
wi aa) - ~ ~ E- 4 Sa 
[*) o om W _ OQ [ae ect 
_ -~ Oo © 4 =. Cdl od ood 
[= ~ = z o aoa xmD 
Lao <= ~- - Ot wy rat ot (S 
uy Y <a = N Fe wl Zz WwW [s Mian ys 4 
wt w= an za it W oa Lal 
be —_ Oo W uo bn x Ire 
ama a mt ULL x x - ~ 0x Enon _ = 
[a] tel aw | ~ mA) > am (5 ad — A Wu Wu ~~ =n 
= - ioe) Am << re = a NS me Ss -o ~t- 3 ~J a OO oO WS 
- 30 ae ia «a IO< -~ §—= Na 7D -~ = NA ~ ud <= a ee 
vw. ge e Zee a 2 ek a x Zz x eet =O am wt 
Sal oO Pa So) WY z “+ Ge S D = ~*~ Oe cad Fee SANS 
oo e ae Om ~ ~ ON >» © es et rN HS OF se wo “et Se) 
— “np OUD DOO a a = (on lent ime # (D> ZzZaee 8D ANN 
= OQ — = cS hH~w OW w i & OF & oO mar NINN & 
ee | Ca “Th Cel oo ae a 2 wm EO Ae ee OS 
ZzwrZz = Zz wZwH= rom O~ eO— YX +O ele OF ce ee CxmINN 
eet oy mae SUS Zeer eZee ns Ze WF He me OZ - mUNOOOAIM & 
Redd & Dried tp Il CRO Ome oe Oe me ee met OO IMO Oh OO ee 
Ww eta ol) em OemL) COW we mete 8 Oe Ae ate sede WE YW TCS xmMx<NL 


LOH W KOUTA NZI oS HOOT tf SUC HE SLU HO HE PLU EN EW Pe WW ah OM Ome 

aod PED Re ZF el ADD ON ME OQDWMeSMS RANMA MDWMSO>Dr 20fp—ee~eTwewe~m 

> '2 $F © £ 9M OXY ZZeHicte ze Zt ze Zr AWS Zee RR eRe Ee 

DOWNY OM IES MS OM Bee 1} EDN RR Rte Ne A eT RE OME REINE Sree sx |) odd 

AAD ASR RZ ZAM HF MOREE ZONZZNZEZAZANZEONZ2NZEZAN2Z6-20 Se SSE TREES TE 
to —2 fiuw OK I~SZzZate Pamet bet en GZ tert Leelee teed ae ee Seed ool en Pe as ee Ae Maa ean A) 
DOKOULDOMALON eM YONOTLE OK LKOYVOULOMOLFOCYOYOLWLCOVCOU ~oOaodCconqg0og002 
CSS ee ee ee 

ont Samad 
Wm Oo OD OOM x on Naw my Sal ane OMDROANMATUY 
Hs « aA ANS N NO N NA NINNN OC MMM 


Oo oO Uo 








~ 
tu 
Q 
2 
= 
De oO 
faa) ~ ' 
be ! 
Ss - 
= « ; 
- ) 
~ ron 
Oo e 
uN a 
= tb 
oO Q 
uw G 
as = 
ton 
= = 
* Oo 
aa x 
mt LEY tL. 
zo t+ 
“oO jes Ge 
Su - 2 
om <r 
zOQ [a - 
muy © ve) 
C~ ~ 
Zoe -_ ~ a i 
Lh WY WY — ' 
a uy _ — ee be 
Som ts Q~ WW 
Own CG = ze, Usk 
Oa j— e -~ » me 
zea ~ _~ mam @, Cet = 
OD boot p Cys . as = it p raee | 
hm x ~ Om e BO el ~ 
Maz es» Gy it << — Za sem 
LU ed we. = ~s 3 mae - eM 
cow us SS om - re | stu 
~au z z wre Cy On rt 
La 4 -* © iu wee a Ze mH a md 
Wu - OoOMnm Zz a1iagw oP a on 
Zoe = On elL as ae az — NON 
ee C st ZORA GB OO wet Ml: & eh ees 
cet ed Dl | ed 20 | VM Mw pas NOAD~~ ~~ 


sev m= BWC it~ ft #2. Ms che pre eZ 

OCS Sie BAO SAHRA DH SH MOOR oy 

KEWZOWMMNAIZS ZFE+ZVMnwse I Ye ZZRESOMED 

DES Tee we Hee OS SO EO 

POOF OOCYSIIYDIYULOOUWLORAODY MONO OWFZ r 
NVOROsM2ZOLYSZEXMZEMZOCHMHOerMY~UGAAOML Waw 


ot 
ot N mat iN OMmO Oo 





INITIAL DISTRIBUTION LIST 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Library 
Naval Postgraduate School 
Monterey, California 93940 


Director, Systems Analysis Division (OP 96) 
Office of the Chief of Naval Operations 
Washington, D. C. 20350 


Department of the Army 
Civil Schools Branch, OPQ, OPD 
Washington, D. C. 20315 


Professor Harold Greenberg 
Department of Operations Analysis 
Naval Postgraduate School 
Monterey, California 93940 


Captain Robert J. Bohls 
346 Ardennes Circle 
Fort Ord, California 93941 


Department of Operations Analysis 


Naval Postgraduate School 
Monterey, California 93940 


46 


No. Copies 


20 





= niente 





UNCLASSIFIED 


Security Classification 


RST SOE Si DI PRR RANT SL 3 TAS AEE PATEL TERESA LS AECL RBS PIE ATT SIRE BG LEY BIE ES UMN ES ELL GSE EEE TERE TORTI LN IE TE TY IIE TL SOLES LL EEE STE IST EES GLE | 


q COCUMENT CONTROL DATA - R&D 


Gocurity clagaification of title, body of abstract and indexing ennotation muet be entered whon the overall raport ta classified) 





RA STENTS EOS 


Za. REPORT SECURITY C LASSIFICA TION 
Unclassified 


SRA PRL AEN RRS IN HD ET PALA IEE BT BEAR LICR TE IL A OLAS OR OR ES SUSE RS SORE 


yp A RIDA II Ti 
pt. “ORIGIHATIN G ACTIVITY (Corporate author) 


Naval Postgraduate School 
Monterey, California 93940 





mera 


oe REPORT TITLE 
On Network Analysis 


| 4. DESCRIPTIVE NOTES (Type of report and inchrelve Gatea) 
Thesis . 

fi cy a 

E 8. AUTHOR(S) (Last name, Hret namo, initial) 

i 4 

j 5 


BOHLS, Robert J., Captain, United States Army 


i 
! 


6. REPORT OATE 


June 1968 = 


acne ret nacnet name 
1 Ga. CONTRACT OR GRANT NO. 








7a. TOTAL NO. OF PAGES 72, NO. OF REFS i 


46 19 


98. ORIGINATOR'S REPORT NUMBER(S) 


+, PROJECT No. 







°. Pinky bo, ; 
9b. Sie mpor ae nat waen app opted io car unt fet eet 
d. its : 
p10. AVAILABILITY/LIMITATION NOTICES ~ pele ; . i 
, if. SUPPLEMENTARY NOTES : 12. SPONSORING MILITARY ACTIVITY } 
; 


Naval Postgraduate School 
i Monterey, California 


13. ABSTRACT 
ie 

} This thesis reviews four of the existing matrix methods for finding 

: the shortest path in a network, including the little known matrix method 

A by Floyd. Floyd's method is then extended to determine all best paths. 

' After a brief review of the nth best path problem, Floyd's method jis 
again extended to determine the nth best path. Finally, the nth best path 
problem is investigated to determine its applicability to the traveling 
salesman problem. 


LRA ATR SES LAT STAI AER TREE SESS SG ASIN. ANIL TI I RUT PAN SEH TSIEN SRDS TOT YY LT LEI TEE 


DD 12". 1473 47 | ROIS Senin ae 


Szeurity Cleasifi cation 


Pak 


_UNCLASSIFIED 


Security Classification 
amar resect 


Shortest Path Prcblem 
All Alternate Be: t Faths Problem 
Nth Best Path Problem 


Traveling Salesman Probiem 





re coe 


GEA =F 88? 


Ee ESR RESTATED EEE LINE LSE TTT 


DD "1473 (Back) 
BPy 


48 


LINK A 


ROLE 





CGE BBLT FEIT PEER. 


UNCLASSIFIED 


Security Classification 


5A cM ADIOS EAB 4 ODS A BET DIE EIT ETT AS LS ESI AO I ES STE 2) SEEAIMENTS TAA A SATE ESTEE 
LINK B LINK ¢€ 


A-31409 





